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Abstract 

The reductive perturbation method has been employed to derive 
the Korteweg-de Vries (KdV) equation for small but finite amplitude 
ion-acoustic waves. The Lagrangian of the time fractional KdV equa- 
tion is used in similar form to the Lagrangian of the regular KdV 
equation. The variation of the functional of this Lagrangian leads to 
the Euler-Lagrange equation that leads to the time fractional KdV 
equation. The Riemann-Liouvulle definition of the fractional deriva- 
tive is used to describe the time fractional operator in the fractional 
KdV equation. The variational-iteration method given by He is used 
to solve the derived time fractional KdV equation. The calculations 
of the solution with initial condition sech(cx)^ are carried out. The 
result of the present investigation may be applicable to some plasma 
environments, such as ionosphere. 
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1 Introduction 

The ion-acoustic solitary wave is one of the fundamental nonlinear wave 
phenomena appearing in fluid dynamics [1] and plasma physics [2, 3]. There 
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are different methods to study nonlinear systems. Washimi and Taniti [4] 
were the first to use reductive perturbation method to study the propagation 
of a slow modulation of a quasi-monochromatic waves through plasma. And 
then the attention has been focused by many authors [5, 6]. The evolution of 
small but finite- amplitude solitary waves, studied by means of the Korteweg- 
deVries (KdV) equation, is of considerable interest in plasma dynamics. 

Because most classical processes observed in the physical world are non- 
conservative, it is important to be able to apply the power of variational 
methods to such cases. A method used a Lagrangian that leads to an Euler- 
Lagrange equation that is, in some sense, equivalent to the desired equation 
of motion. Hamilton's equations are derived from the Lagrangian and are 
equivalent to the Euler-Lagrange equation. If a Lagrangian is constructed 
using noninteger-order derivatives, then the resulting equation of motion can 
be nonconservative. It was shown that such fractional derivatives in the 
Lagrangian describe nonconservative forces [7, 8] . Further study of the frac- 
tional Euler-Lagrange can be found in the work of Agrawal [9, 10], Baleanu 
and coworkers [11, 12] and Tarasov and Zaslavsky [13, 14]. During the last 
decades. Fractional Calculus has been applied to almost every field of science, 
engineering and mathematics. Some of the areas where Fractional Calculus 
has been applied include viscoelasticity and rheology, electrical engineering, 
electrochemistry, biology, biophysics and bioengineering, signal and image 
processing, mechanics, mechatronics, physics, and control theory [15]. 

To the author's knowledge, the problem of time fractional KdV equation 
in coUisionless plasma has not been addressed in the literature before. So, 
our motive here is to study the effects of time fractional parameter on the 
electrostatic structures for a system of coUisionless plasma consisting of a 
mixture of warm ion-fiuid and isothermal electrons. We expect that the 
inclusion of time fractional parameter will change the properties as well as 
the regime of existence of solitons. 

Several methods have been used to solve fractional differential equations 
such as: the Laplace transform method, the Fourier transform method, the 
iteration method and the operational method [16]. Recently, there are some 
papers deal with the existence and multiplicity of solution of nonlinear frac- 
tional differential equation by the use of techniques of nonlinear analysis 
[17-18]. In this paper, the resultant fractional KdV equation will be solved 
using a variational-iteration method (VIM) firstly used by He [19]. 

This paper is organized as follows: Section 2 is devoted to describe the for- 
mulation of the time-fractional KdV (FKdV) equation using the variational 
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Euler-Lagrange method. In section 3, variational-Iteration Method 9 VIM) is 
discussed. The rcsuhant time-FKdV equation is solved approximately using 
VIM. Section 5 contains the results of calculations and discussion of these 
results. 

2 Basic Equations and KdV Equation 

Consider a coUisionless ionization-free unmagnetized plasma consisting of a 
mixture of warm ion-fluid and isothermal electrons. The basic equations 
describing the system in dimensionless variables is governed by [6] : 

0, (l.a) 
0, (l.b) 

0, (l.c) 

0, (l.d) 
0. (l.e) 

In the earlier equations n{x,t) and ne{x,t) are the densities of ions and 
electrons respectively, u{x,t) is the ion flow velocity, p{x,t) is the ion pres- 
sure, t) is the electric potential, x is the space co-ordinate and t is the 
time variable, a = Ti/Tf. « 1 is the ratio of the ion temperature Tj to 
the electron temperature Tg. All these quantities are dimensionless, being 
normalized in terms of the following characteristic quantities: n{x, t) and 
ne{x, t) by the unperturbed electron density Uq, u{x, t) by the sound velocity 
{KTe/miY^'^ ; p{x,t) and (f){x,t) by noKTi and KT^je, respectively. The 
time variable t and the space x by the inverse of the plasma frequency ijJ~1 — 
{Ane'^no /■mi)~^^'^ and the electron Debye length Xd = {KT^/ Ane^noY^'^ , re- 
spectively. K is the Boltzmann's constant and rrii is the mass of plasma 
ion. 

According to the general method of reductive perturbation theory, we 
introduce the stretched variables 

T = e-H, ^ = e-^{x-vt), (2) 
3 



dn{x,t) d[n{x,t)u{x,t)] 



dt 
a 



+ 



dx 



dt 



dx 



n{x, t) dx 
-}p{x,i 
d'^(t){x,t) 



d d 



dp{x,t) d(j){x,t) 
dx 
du{x, t) 



dx' 



+ n{x,t) — ne{x,t) — 
ne{x,t) - exp[(f){x,t)] = 



where V is the unknown phase velocity. All the physical quantities appeared 
in (1) are expanded as power series in e about the equilibrium values as 



n{C,T) = l + eni(e,T) + eV(e,T)+ (3a) 

u{C,t) = eiii(e,T) + eV(e,T)+ (3b) 

p{C,t) = l + epi(e,T)+e'p2(e,T)+ (3c) 

<P{C,r) = eM^,r) + e'MC,r)+ (3d) 

where e represents the amplitude of the perturbation. We impose the bound- 
ary conditions that as |^| oo, n — Ue = p = 1, u = (f) — 0. Substituting 
(2) and (3) into the system of equations (1) and equating coefficients of hke 
powers of e. Then, from the lowest, second-order equations in e and by elim- 
ination the second order perturbed quantities n2{^,T), U2{^,t), P2{^,t) and 
02(C)''")) '^^ obtain the following KdV equation for the first-order perturbed 
potential: 



^MC, r)+A M^, t)-^M^, r) + B ^^C, r) = 0, (4.a) 

where both the nonlinear coefficient A and the dispersion coefficient B are 
given by 



In equation (4. a), 0i(.C, r) is a field variable, ^ is a space coordinate in the 
propagation direction of the field and r e T (= [0,To]) is the time variable. 
The resultant KdV equation (4. a) can be converted into time-fractional KdV 
equation as follows: 

Using a potential function V{^^ r) where <f>i(C, r) = V^(C, r) gives the 
potential equation of the regular KdV equation (1) in the form 



V^,{i, t) + A V^{i, r)^55(e, t) + B %^^(e, r) = 0, (5) 

where the subscripts denote the partial differentiation of the function with 
respect to the parameter. The Lagrangian of this regular KdV equation (4. a) 
can be defined using the semi- inverse method [20, 21] as follows. 

The functional of the potential equation (5) can be represented by 
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J{V) = J J dT{V{^,T)[c^V^r{^,T)+C2AV^{^,TX^{^,T) 
R T 

+C3i?%«(C,r)]}, (6) 

where ci, C2 and C3 are constants to be determined. Integrating by parts and 
taking Vt\ji = V^Ir — V^\t = lead to 



J{V) = Jd^jdT{V{C,T)[-C^V^{C,T)Vr{C,T)-^C2AVl{C,T) 
R T 

+C3W4(e,r)]}. (7) 

The unknown constants Cj {i = 1,2,3) can be determined by taking the 
variation of the functional (7) to make it optimal. Taking the variation of this 
functional, integrating each term by parts and make the variation optimum 
give the following relation 

2ciV^^{^, r) + 3c2AV^{^, r)%(e, r) + 2csBV^^^^{^, r) = 0. (8) 

As this equation must be equal to equation (5), the unknown constants 
are given as 

ci = 1/2, C2 = 1/3 and C3 = 1/2. (9) 

Therefore, the functional given by (7) gives the Lagrangian of the regular 
KdV equation as 



m, = --V^{C,t)VAC,t) - -AVl{i,T) + -W|fer). (10) 

Similar to this form, the Lagrangian of the time-fractional version of the 
KdV equation can be written in the form 



FioD^V, V^, %) 



1 



< 



2 

a < 1, 



(11) 
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where the fractional derivative is represented, using the left Riemann-Liouville 
fractional derivative definition as [16] 



k-1 < a<l,t e[a,b]. (12) 
The functional of the time-FKdV equation can be represented in the form 

J{V) - Jd^J drFioD^V, 1^,%), (13) 

R T 

where the time-fractional Lagrangian F{qD"V, V^, V^^) is defined by (11). 

Following Agrawal's method [9, 10], the variation of functional (13) with 
respect to V{^,t) leads to 

SJiV) M^WV + §-^6V, + §-W,,}. (14) 

R T 

The formula for fractional integration by parts reads [9, 16] 



f dtf{t)aD^g{t) = f dtg{t)tDtf{t), f{t), g{t) e [a, b], (15) 

J a J a 

where t^*^, the right Riemann-Liouville fractional derivative, is defined by 
[16] 



k-1 < a<l,t e[a,b]. (16) 
Integrating the right-hand side of (14) by parts using formula (15) leads 



to 



R T 
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where it is assumed that SV\t = SV\ji = SV^Ir = 0. 

Optimizing this variation of the functional J{V), i. c; 6J{V) = 0, gives 
the Euler-Lagrange equation for the time-FKdV equation in the form 



Substituting the Lagrangian of the time-FKdV equation (11) into this 
Euler-Lagrange formula (18) gives 

-^ri^Tom^) + ^oi^."T^e(e,r)+>iye(e,r)%(e,r)+W^^^^(e,r) = 0. (19) 

Substituting for the potential function, V^(C,r) = 0;^(^, r) = $(^, r), 
gives the time-FKdV equation for the state function $(^, r) in the form 

^ [oD^m t) r)] + A^{i, r)^^{i, r) + ^^^(e, r) = 0, (20) 

where the fractional derivatives o-D" and tD^^ are, respectively the left and 
right Riemann-Liouville fractional derivatives and are defined by (12) and 
(16). 

The time-FKdV equation represented in (20) can be rewritten by the 
formula 

\ ^D^$(e, t) + A *(e, r)$^(e, t) + B %^(e, r) = 0, (21) 

where the fractional operator ^D" is called Riesz fractional derivative and 
can be represented by [10, 16] 



d 



[f dT\t-r\'—'f{r)], 

J a 



2 r(A; - a) dt^ ' 

k-l < a<l,te[a,h]. (22) 

The nonhnear fractional differential equations have been solved using dif- 
ferent techniques [16-20] . In this paper, a variational-iteration method (VIM) 
[21, 22] has been used to solve the time-FKdV equation that formulated using 
Euler-Lagrange variational technique. 
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3 Variational-Iteration Method 



Variational- iteration method (VIM) [21] has been used successfully to solve 
different types of integer nonhnear differential equations [22, 23]. Also, VIM 
is used to solve hnear and nonlinear fractional differential equations [24, 25]. 
This VIM has been used in this paper to solve the formulated time-FKdV 
equation. 

A general Lagrange multiplier method is constructed to solve non-linear 
problems, which was first proposed to solve problems in quantum mechanics 
[21]. The VIM is a modification of this Lagrange multipher method [22]. The 
basic features of the VIM arc as follows. The solution of a linear mathemat- 
ical problem or the initial (boundary) condition of the nonlinear problem 
is used as initial approximation or trail function. A more highly precise 
approximation can be obtained using iteration correction functional. 

Considering a nonlinear partial differential equation consists of a linear 

part LU {x, t), nonlinear part NU {x, t) and a free term f{x, t) represented as 

LU(x,t) + NU(x,t) ^ f(x,t), (23) 

where L is the linear operator and N is the nonlinear operator. According 
to the VIM, the (n+ l)th approximation solution of (23) can be given by the 
iteration correction functional as [21, 22] 



Un+l{x,t) = Un{x,t)+ [ dTX{T)[LUn{x,T)+NUn{x,T)-f{x,T)], 

Jo 

n > 0, (24) 

where A(r) is a Lagrangian multiplier and tJn{x,T) is considered as a re- 
stricted variation function, i. e; 5?7„(a;,r) = 0. Extreme the variation of 
the correction functional (24) leads to the Lagrangian multiplier A(r). The 
initial iteration can be used as the solution of the linear part of (23) or the 
initial value U{x, 0). As n tends to infinity, the iteration leads to the exact 
solution of (23), i. e; 

U{x,t) = lim Un{x,t). (25) 

n— >-oo 

For linear problems, the exact solution can be given using this method 
in only one step where its Lagrangian multiplier can be exactly identified. 
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4 Solution of the time-FKdV equation 



The time-FKdV equation represented by (21) can be solved using the VIM 
by the iteration correction functional (24) as follows: 

Affecting from left by the fractional operator on (21) leads to 



9t ^"'^ u-r -v-'^i--r(a-l) 

- ^D^'--[A r)-$(e, r) + B r)], 
< a<l, re[0,To], (26) 
where the following fractional derivative property is used [16] 

k-l < P<k. (27) 
Riesz 

fractional integral ^/^~° that is defined by [10, 16] 



As a < 1, the Riesz fractional derivative qD^°' ^ is considered as Riesz 



Olr" fit) = l^rVit) + rhVm 



1 1 



2r(a) 



I cir|t-r|"-V(r), a > 0, (28) 

J a 



where oI^°'f{t) and ^/^"/(t) are the left and right Riemann-LiouvuUe frac- 
tional integrals, respectively [16]. 

The iterative correction functional of equation (26) is given as 



$„+l(e,T) = $„(^,r) + ^'c^T'A(T'){^$n(e,T') 



-o''/.V"*n(e,r')|r'=0 



r{a-i) 

+ X^"[A $.(e, r )^$n(e, r') + i?^$n(e, r')]},(29) 
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where n > and the function $„(^, r) is considered as a restricted variation 

function, i. e; 5$„(^,r) = 0. The extreme of the variation of (29) using the 
restricted variation function leads to 



dr'Xir') d—<^n{C,r') 

= S^niL r) + A(r) 6^n{L r) - J^^ dr'-^Xir') r') = 0. 

This relation leads to the stationary conditions 1+A(r) — and ^A(t') = 
0, which leads to the Lagrangian multiplier as A(t') = — 1. Therefore, the 
correction functional (29) is given by the form 



$„+i(e,r) = $„(e,r)-^"dT'{^$„(e,r') 



o'^/.r"$„(e,r')| 



la-2 



T(a- 1) 
d 



+ SDr)-'[A *„R, r')^4>„(«, r') + B ^*„(6 t')1}(30) 

where n > 0. 

In Physics, if r denotes the time-variable, the right Riemann-Liouville 
fractional derivative is interpreted as a future state of the process. For this 
reason, the right-derivative is usually neglected in applications, when the 
present state of the process does not depend on the results of the future 
development [9]. Therefore, the right-derivative is used equal to zero in the 
following calculations. 

The zero order correction of the solution can be taken as the initial value 
of the state variable, which is taken in this case as 

^o(e,T) =$(e,0) = Asech2(cO. (31) 

where = ^ and c = |-\/^ constants. 

Substituting this zero order approximation into (30) and using the defi- 
nition of the fractional derivative (22) lead to the first order approximation 
as 
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$i(e,r) = Aosech\cO 

+2AqC sinh(cO sech^(c^) 

*[Ac'B + {AoA - Uc'B) sec /^^(cO] , ■ (32) 

Substituting this equation into (30), using the definition (22) and the 
Maple package lead to the second order approximation in the form 



$2(e,r) = Aosech2(cO 

+2AoC sinh(c^) sech(c^)^ 

*[4c^B + {AqA - 12c^B) sech(cO^ 



+2AoC^ sech(c,^)2 

*[32c'^B^ + 16c^B(5AoA - 63c^B) sech(cO^ 
+2{3AlA^ - 176Aoc'^AB + IGSOc^S^) sech(cO^ 

-TiAlA^ - A2Aoc^AB + 360c^5^) sech(c{)^]— ^ 



2a 



T{2a + 1) 

+4yloC^ sinh(c^) scch(cO^ 
*[32c^B^ + 2Ac^B{AoA - Uc^B) sech(cO^ 
+A{AlA^ - 32Aoc^AB + 2A0c^B^) sech(ce)^ 
-5{AlA^ - 2AAqc^AB + lAAc^B^) sech(cO^] 

T{2a + 1) r^" 
*[r(a + l)]2r(3a + l) 



(33) 



The higher order approximations can be calculated using the Maple or 
the Mathematica package to the appropriate order where the infinite approx- 
imation leads to the exact solution. 



5 Results and calculations 

For a small amplitude ion-acoustic solitary wave in unmagnetized coUisionless 
plasma consisting of a mixture of warm ion-fluid and isothermal electrons, 
we have obtained the Korteweg-de Vries equation by using the reductive 
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perturbation method. The time fractional Korteweg-de Vries equation has 
been derived using the variational technique [9, 20]. The Riemann-Liouvulle 
fractional derivative is used to describe the time fractional operator in the 
time-FKdV equation. He's variational-iteration method [21, 22] is used to 
solve the derived time-FKdV equation. The calculations in this work is 
carried out using the fourth order approximation of the VIM. 

However, since one of our motivations was to study effects of the ratio of 
the ion temperature to the electron temperature o and the time fractional 
order a on the existence of solitary waves. The present system admits a 
solitary wave solution for any order approximation. In Fig (1), a profile 
of the solution of the time FKdV equation as a sohtary pulse is obtained. 
Figure (2) shows that the both the amplitude and the width of the solitary 
wave decrease with the increase of the temperatures ratio a . Also, the time 
fractional order a and the velocity v increase the soliton amplitude as shown 
in Fig (3). In summery, it has been found that amplitude of the ion-acoustic 
waves as well as parametric regime where the solitons can exist is sensitive 
to the temperatures ratio a and the time fractional order a. Moreover, the 
time fractional order a plays a role of higher order perturbation theory in 
increasing the soliton amplitude. The application of our model might be 
particularly interesting in some plasma environments, such as ionosphere. 
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6 Figure Captions 

Fig. (1): The electrostatic potential $(^, r) against the position ^ and the 
time T ior a — 0.5, v — 0.4 and a — 0.5. 

Fig. (2): The electrostatic potential $(^, r) against the position ^ and 

the temperatures ratio a at r = 1, ?; = 0.4 and a = 0.5. 

Fig. (3): The electrostatic potential $(^, r) against the velocity v and 
the fractional order a at ^ = 0, t = 1 and a = 0.5. 
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